
## From TODO:

m01 <- matrix(data=c(11.0, 12.0, 13.0, 14.0, 15.0,
                21.0, 22.0, 23.0, 24.0, 25.0,
                31.0, 32.0, 33.0, 34.0, 35.0,
                41.0, 42.0, 43.0, 44.0, 45.0,
                51.0, 52.0, 53.0, 54.0, 55.0,
                61.0, 62.0, 63.0, 64.0, 65.0),
              nrow=6, ncol=5,byrow=TRUE)

m01
mtm <- t(m01) %*% m01
mtm

qr(mtm)

solve(mtm)

